Finite-size scaling of correlation functions in 
one-dimensional Anderson-Hubbard model 



0^ 
O 
O 
(N 

D 



Satoshi Nishimoto 

Leibmz-Institut fiir Festkorper- und Werkstoffforschung Dresden, D-01171 Dresden, Germany 

Tomonori Shirakawa 
Department of Physics, Chiba University, Chiba 263-8522, Japan 
Institut fiir Theoretische Physik, Leibniz Universitdt Hannover, D-30167 Hannover, Germany 

(Dated: February 19, 2009) 

We study the one-dimensional Anderson-Hubbard model using the density-matrix renormalization 
group method. The influence of disorder on the Tomonaga-Luttinger-liquid behavior is quantita- 
tively discussed. Based on the finite-size scaling analysis of density-density correlation functions, we 
find the following results: i) the charge exponent is significantly reduced by disorder at low filling 
and near half filling, ii) the localization length decays as ^ ~ A~^, where A is the disorder strength, 
independently of the on-site Coulomb interaction as well as band filling, and iii) the localization 
length is strongly suppressed by the on-site Coulomb interaction near half filling in association with 
the formation of the Mott plateaus. 

PACS numbers: 71.10.Fd, 71.10.Pm, 71.23.-k, 71.15.Dx 
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We know that lattice disorder or defects in the elec- 
tron system lead to various intriguing effects. One of 
the most famous examples is Anderson localization [l[ ; a 
large number of phenomena is explained in terms of the 
localization of quantum states Meanwhile, it is widely 
recognized that a full understanding of the experimen- 
tally observed phenomena has a need for taking correla- 
tion effects, as well as disorder, into account. Hence, the 
Anderson localization in the presence of electron-electron 
interactions has been an inevitable problem in the field of 
condensed matter physics [3||. A hot topic is the appear- 
ance of "zero bias anomaly [J]", i.e., a suppression of the 
density of states A{E) at the Fermi energy Ef, by the 
interplay between disorder and short-range Coulomb in- 
teraction 0, 0] . Recently, it was reported that the similar 
anomaly occurs even in one-dimensional (ID) system Q. 
It is possibly related to a vanishing of the photoemission 
spectral weight around Ep (or "pseudogap-like" behav- 
ior) observed in quasi- ID charge-transfer materials TTF- 
TCNQ and Bechgaad salts (TMTSF)2X 

Generally, ID systems play a key role for elucidating 
the physical properties of solids because many features 
of the electronic states can be discussed rigorously. It 
has been also confirmed that the ground state and low- 
lying excitations can be described within the framework 
of the Tomonaga-Luttinger-liquid (TLL) theory [l^] even 
in the presence of both disorder and electron-electron in- 
teraction 13, U I ■ Therefore, we are allowed to measure 
the effect of disorder as a modification of the TLL states. 
The properties of TLL are characterized by a few quan- 
tities; most notably, the charge exponent Kp determines 
the long-range behavior of various correlation functions 
in the ID metallic state. In this regard, a point to ponder 
is that the all eigenstates of a disordered ID system are 
exponentially localized in the asymptotic sense • 



In this Letter, we study a ID Hubbard model in 
the presence of disorder, namely ID Anderson-Hubbard 
(AH) model. Using the density-matrix renormalization 
group (DMRG) technique the system-size depen- 

dent density-density correlation functions are calculated 
for various on-site Coulomb interactions and disorder 
strengths. We then propose a finite-size scaling method 
of the correlation functions for obtaining the modified 
TLL charge exponent and localization length of the elec- 
trons. In consequence, a quadratic decay of the localiza- 
tion length with the inverse of disorder strength is con- 
firmed for any interaction and band filling. We also find 
a drastic reduction of the charge exponent with disorder 
except around quarter filling and a strong suppression of 
the localization length with the on-site Coulomb interac- 
tion near half filling. 

The Hamiltonian of the ID AH model is written as 



L-l 



i=l 



(1) 

where Cia is annihilation operator of an electron with 
spin a {—^ or J,) at site i, L is the system length, and 
nia = cl^Cia is number operator. The nearest-neighbor 
hopping integral t and on-site Coulomb interaction U are 
assumed to be constant over the system. The random 
on-site potential Si {i ~ 1,---,L) is defined by a box 
probability distribution V{ei) ~ 0{l/2 — \ei\) with the 
step function 9(x) and the disorder strength is controlled 
by A. The band filling is n = N/L where TV is the 
total number of electrons. We set t = 1 as energy unit 
hereafter. 

We consider the long-range behavior of the density- 
density correlation in the presence of disorder. We focus 
on the case of C/ > and < A < C//2, where the sys- 
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tern is in the localized phase. Thus, the density-density 
correlation functions may be defined like [11, [3 



C{r) 



= e 



C{r) 



A=0 



(2) 



with the asymptotic behavior in the absence of disorder 



C(r)|A=o = - 



Acos(2fcFr) _3/2 

In ' (r)H , (3) 



rl + Kp 



where ^ corresponds to the localization length of the elec- 
trons. Calculating the Fourier transformation of Eq.Q, 



C(<7i) = ^C(r 



(4) 



with qi = 27r/L, we obtain 



K*„ - 1 



27r2 ^ 



1 



(5) 



for gi <C 1 (L ^ 1). The parameter K* is interpreted as 
a modified TLL charge exponent due to the disorder. In 



the limit of weak disorder, 



is equivalent to Kp and ^ 
17| . we cal- 



tends to be infinite. By the DMRG method 
culate C{qi) on random sampling 300 (500) realizations 
of ^(eO for L = 128, 112, 96, 80, 64, 48, (32, and 16); 
then, we take an average of the results for each system 
size in order to obtain physically meaningful values of 
C{qi) [The averaged value is denoted as C{qi)]. For ac- 
curate calculation, the open-boundary conditions (OBC) 
are applied and we keep up to m = 2400 density-matrix 
eigenstates in the renormalization procedure. Note that 
special attention should be paid to the convergence of the 
calculation because the DMRG wave function is apt to 
get trapped in a 'false' ground state for disordered sys- 
tem. Thus, K* and ^ will be estimated by fitting our 
numerical data of C{qi) with Eq.([n]). 

For instance, we demonstrate the fitting of C{qi) at 
U = 10 and n = 1/2. Figure[lja) shows the averaged val- 
ues of C{qi) as a function of 1/(L — (5) for various disorder 
strengths A. Since a certain quantity of the correlations 
is missing at system edges due to the OBC, L should be 
replaced with L — 6 {5 > 0) for an excellent finite-size- 
scaling analysis. The correction factor S, which would be 
of the order of the lattice spacing, is determined to re- 
produce the relation K* = Kp in the absence of disorder; 
as expected, we obtain a quasi-infinite localization length 
(l/l^l < 1.0 X 10-5) with 5 « 0.8 - 1.2 at A = for aU 
fillings and interaction strengths considered in this study. 
In this way, the DMRG results of C{qi) are well fitted by 
Eq.([5]) even for finite A values, as shown in Fig.[T](a). The 
estimated values of K* and ^ are plotted as a function of 
A in Fig.[T](b) and (c), respectively. The exponent K* is 
reduced with increasing A because the forward scatter- 
ing processes are enhanced by the disorder [lij . However, 



(a) 0.08 

0.06 
5 0.04 
0.02 



— ' 1 ' — 

A = 

1 

2 
3 

,r-';>.--''.'-' 
»*■''.■'.■■''--' 



0.02 0.04 0.06 0.08 




FIG. 1: (a) Finite-size-scaling analyses of the charge cor- 
relation function C'(gi) for various disorder strengths A at 
U = 10. The dotted lines are the fitting with Eq. ((5}. (b) 
Fitting results of the TLL charge exponent K*p as a function 
of A. The dashed line is a second-order polynomial in A. (c) 
Fitting results of the localization length ^ as a function of A. 
Insets: Log- log plot of ^ with A (open circles). The values of 
^ estimated from the Drude weight are also plotted as solid 
circles. 



the shortening of the localization length ^ on A seems to 
be much more drastic. Wc find that ^ decays as a power 
law with A, i.e., i = CoA"", if A < C//2. At n 1/2, we 
estimate ^o ~ 950 and a w 2 [sec the inset of Fig. [T](c)]. 

In order to check the accuracy of the above analysis, 
we also estimate the localization length from an expo- 
nential decay of the Drude weight D(V) with the system 
length L [l^i Ell and perform a cross-check. The Drude 
weight is calculated on random sampling 100 realizations 
of Vie.i) for 64, 48, 32, and 16 using a recently proposed 
method [20| . The periodic boundary conditions are ap- 
plied for this calculation. By fitting of the averaged re- 
sults of D{V) [= D{V)\ with a formula 



5(L)-7exp(--), 



7 > 0, (6) 
we can obtain ^ without any difficulty. As seen in the in- 
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FIG. 2: Modified charge exponent K*p as a function of the 
band filhng n for various values of the disorder strength A at 
(a) J7 = 10 and (b) U = 2. The dashed hues are guides to the 
eye. 



FIG. 3: The locahzation length ^ as a function of the disorder 
strength A for (a) ?7 = 10 and (b) [7 = 2. The dashed lines 
are the fitting function ^ = for small A. 



set of Fig. [Ijc) , there is an excellent agreement between 
the resuhs from D{L) and C{qi). Therefore, we can con- 
firm the validity of our method. 

Let us now see how the charge exponent is modified 
by the disorder. Figure [2] shows the values of K* ob- 
tained by the above fitting procedure as a function of 
the band filling n for various disorder strengths A at 
(a) strong {U = 10) and (b) weak (C7 = 2) interaction 
strengths. In the case oi U = 10, we see that the charge 
exponent is most severely affected by the disorder near 
n = 0; namely, K* is drastically reduced by A since 
most of (or all) the electrons are easily trapped at sites 
having low on-site potential. Also near n = 1, K* de- 
creases rapidly with increasing A. This can be explained 
in terms of the formation of Mott plateaus assisted by the 
disorder 2]|. In other words, the system will be essen- 
tially in the Mott insulating state for tiny doping since 
the holes are strongly localized. On the other hand, the 
effect of disorder appears to be relatively weak around 
quarter filling because significant "untrapped" particles 
still remain and they could transfer. Consequently, the 
charge exponent is hardly affected for small A 1) at 
n ~ 0.3 — 0.7. It is consistent with a rather unstable Mott 
plateau far away from half filing 2l|. We then turn to 



from half filing and, however, the effect of A near half 
filling seems to be weaker in comparison with the case of 
U = 10. It would be concerned with the fact that the 
formation of the Mott plateau is harder for smaller U. 

Of particular interest is the evolution of the localiza- 
tion length ^ with increasing the disorder strength A. 
Figure [3] shows the fitting values of ^ as a function of 
A for several kinds of the band filling n at (a) U = 10 
and (b) U = 2. We find that the localization length al- 
ways decay quadratically with the inverse of the disorder 
strength, i.e.. 



A2 



(7) 



the case ofU = 2. Qualitatively the same behavior of K* 
with A is observed as long as the band filling is far away 



when the disorder is weak. The same relation has 
been already proposed for spinless fermions case {U = 
oo) [13, HI], which is the same model as what proposed 
by Anderson [ij. Therefore, the universality of Eq.([7]) is 
confirmed in a weakly disordered system independently 
of the on-site Coulomb interaction and band filling. We 
note that the localization length begins to deviate from 
Eq.(l7]) at relatively small disorder A w 1.5 — 2 near 
n = 0,1 for U = 10. It is because a part of the sites 
is doubly occupied under the strong disorder and the 
density-density correlation no longer obeys Eq.([2]). 

A qualitative determination of the decay ratio in 
Eq.([7]) is another challenging problem. So far, it has been 
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a = {Kp + — 2)/4 is still valid in the presence of the 
disorder, the strong reduction of Kp would be compatible 
with the pseudogap-like behvbior. 
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FIG. 4: Fitting values of the decay ratio for various inter- 
action strengths. 



studied perturbatively [2^ and numerically 1^, 24, 25 1 
only for spinless fcrmions. Figure 0] shows the fitting val- 
ues of ^0 as a function of the band filling n at C/ = 0, 2, 
10, and oo. When the interaction strength is varied from 
U = to 2, ^0 increases for all fillings. It seems to be 
consistent with a prediction that the persistent currents 
are enhanced by the repulsive interactions [2^ and, how- 
ever, ^0 decreases with further increasing U . Especially 
near n = 1, the reduction of with U is quite rapid in 
connection to the stabilization of the Mott plateaus. On 
the other hand, surprisingly, is almost unchanged at 
n < 0.5 for [/ > 2. Consequently, becomes symmetric 
about 11 ~ 1/2 at large C/, which reflects the particlc-holc 
symmetry in the limit oi U — s- cxd. 

In summary, we consider the ID Anderson-Hubbard 
model using the DMRG technique. A finite-size scal- 
ing method of the density-density correlation functions 
for obtaining the modified TLL charge exponent and lo- 
calization length is demonstrated. As results, we find a 
quadratic decay of the localization length with the in- 
verse of disorder strength for any interaction and band 
filling, a drastic reduction of the TLL charge exponent 
with disorder except around quarter filling and, a strong 
suppression of the localization length with the on-site 
Coulomb interaction near half filling. 

Lastly, we make a short remark regarding an explana- 
tion of the pseudogap-like behavior and very small TLL 
exponent Kp « 0.23 ^ observed in the TMTSF salts, 
which is quarter-filled system with dimerization. Re- 
garding the dimer (TMTSF)2 as a site, it may be re- 
duced to a half-filled system. Thus, we will expect a 
strong reduction of the TLL exponent even by tiny dis- 
order. Moreover, if the relation A{E) ^ \E — Ey\°' with 



[1] P.W. Anderson, Phys. Rev. 109, 1492 (1958). 

[2] For a review, B. Kramer and A. MacKinnon, Rep. Prog. 

Phys. 56, 1469 (1993). 
[3] For a review, M. Lewenstein, A. Sanpera, V. Ahufin- 

ger, B. Damski, A. Sen, and U. Sen, Adv. Phys. 56, 243 

(2007). 

[4] A.L. Efros and B.I. Shklovskii, J. Phys. C 8, L49 (1975). 
[5] S. Chiesa, P.. Chakraborty, W.E. Pickett, and 

R.T. Scalettar. larXiv:0804 .4463vl. 
[6] Y. Song, S. Bulut, R. Wortis, W.A. Atkinson, 

|arX iv:08m3_356 Vl. 
[7] H. Shinaoka and M. Imada, Phys. Rev. Lett. 102, 016404 

(2009). 

[8] V. Vescolil, F. Zwick, W. Henderson, L. Degiorgi, 

M. Grioni, G. Gruner, and L.K. Montgomery, Eur. Phys. 

J. B 13, 503 (2000). 
[9] M. Sing, U. Schwingenschlogi, R. Claessen, M. Dressel, 

and G.S. Jacobsen, Phys. Rev. B 67, 125402 (2003). 
[10] D. Loss, Phys. Rev. Lett. 69, 343 (1992). 
[11] A.O. Gogolin and N.V. NV Prokof'ev, Phys. Rev. B 50, 

4921 (1994). 
[12] J. Solyom, Adv. Phys. 28, 201 (1979). 
[13] E. Abrahams, P.W. Anderson, D.C. Licciardello, and 

T.V. Ramakrishnan, Phys. Rev. Lett. 42, 673 (1979). 
[14] S.R. White, Phys. Rev. Lett. 69, 2863 (1992); Phys. Rev. 

B 48, 10345 (1993). 
[15] T. Saso, Y. Suzumura, and H. Fukuyama, Prog. Theor. 

Phys. Suppl. 84, 269 (1985). 
[16] T. Giamarchi and H.J. Schuiz, Phys. Rev. B 37, 325 

(1988). 

[17] S. Ejima, F. Gebhard, and S. Nishimoto, Europhys. Lett. 
70, 492 (2005). 

[18] C.L. Kane and M.P.A. Fisher, Phys. Rev. Lett. 68, 1220 
(1992). 

[19] P. Schmitteckert, T. Schulze, C. Schuster, P. Schwab, and 
U. Eckern, Phys. Rev. Lett. 80, 560 (1998). 

[20] T. Shirakawa and E. Jeckelmann, in preparation. 

[21] M. Okumura, S. Yamada, N. Taniguchi, and M. Machida, 
Phys. Rev. Lett. 101, 016407 (2008). 

[22] D.J. Thouless, J. Phys. C 5. 77 (1972). 

[23] D.C. Herbert and R. Jones, J. Phys. C 4, 1145 (1971). 

[24] G. Czycholl, B. Kramer, and A. MacKinnon, Z. Phys. B 
43, 5 (1981). 

[25] M. Kappus and F.J. Wegner, Z. Phys. B 45, 15 (1981). 
[26] T. Giamarchi and B.S. Shastry, Phys. Rev. B 51, 10915 
(1995). 

[27] M. Dressel, A. Schwartz, G. Griiner, and L. Degiorgi, 
Phys. Rev. Lett. 77, 398 (1996). 



